Algebraic interface-based coarsening AMG preconditioner for multi-scale sparse matrices with applications to radiation hydrodynamics computation

نویسندگان

  • Xiaowen Xu
  • Zeyao Mo
چکیده

Coarsening is one of the crucial components of an algebraic multigrid (AMG) method for iteratively solving sparse linear systems arising from scientific and engineering applications. It largely determines the complexity of the AMG iteration operator. Usually, high operator complexities lead to fast converge of the AMG method, however, they require additional memory and as such do not scale as well in parallel computation. On the other hand, though low operator complexities improve parallel scalability, they often lead to deterioration in convergence. This paper introduces a new type of coarsening strategy, called algebraic interface based coarsening that yields better balance between convergence and complexity for a class of multi-scale sparse matrices. Numerical results for various model-type problems and a radiation hydrodynamics practical application are provided to show the effectiveness of the proposed AMG solver.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Reducing Complexity in Parallel Algebraic Multigrid Preconditioners

Algebraic multigrid (AMG) is a very efficient iterative solver and preconditioner for large unstructured linear systems. Traditional coarsening schemes for AMG can, however, lead to computational complexity growth as problem size increases, resulting in increased memory use and execution time, and diminished scalability. Two new parallel AMG coarsening schemes are proposed, that are based on so...

متن کامل

A fast algebraic multigrid preconditioned conjugate gradient solver

This work presents a new approach for selecting the coarse grids allowing a faster algebraic multigrid (AMG) preconditioned conjugate gradient solver. This approach is based on an appropriate choice of the parameter a considering the matrix density during the coarsening process which implies in a significant reduction in the matrix dimension at all AMG levels. 2005 Elsevier Inc. All rights rese...

متن کامل

On the Scalability of the Albany/FELIX first-order Stokes Approximation ice Sheet Solver for Large-Scale Simulations of the Greenland and Antarctic ice Sheets

We examine the scalability of the recently developed Albany/FELIX finite-element based code for the first-order Stokes momentum balance equations for ice flow. We focus our analysis on the performance of two possible preconditioners for the iterative solution of the sparse linear systems that arise from the discretization of the governing equations: (1) a preconditioner based on the incomplete ...

متن کامل

Fast hierarchical solvers for sparse matrices

Sparse linear systems are ubiquitous in various scientific computing applications. Inversion of sparse matrices with standard direct solve schemes are prohibitive for large systems due to their quadratic/cubic complexity. Iterative solvers, on the other hand, demonstrate better scalability. However, they suffer from poor convergence rates when used without a preconditioner. There are many preco...

متن کامل

DELFT UNIVERSITY OF TECHNOLOGY REPORT 13-11 Evaluation of Multilevel Sequentially Semiseparable Preconditioners on CFD Benchmark Problems Using IFISS

This paper studies a new preconditioning technique for sparse systems arising from discretized partial differential equations (PDEs) in computational fluid dynamics (CFD), which exploit the multilevel sequentially semiseparable (MSSS) structure of the system matrix. MSSS matrix computations give a data-sparse way to approximate the LU factorization of a sparse matrix from discretized PDEs in li...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:
  • Numerical Lin. Alg. with Applic.

دوره 24  شماره 

صفحات  -

تاریخ انتشار 2017